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Measurements  of  Antenna  Pulse  Response 


1.  INTRODUCTION 

The  objective  of  the  work  described  in  the  paper  is  to  obtain  the  pulse  response 
of  antennas.  This  is  accomplished  by  measuring  the  steady-state  response  as  a 
function  of  frequency,  transforming  to  low  frequency,  multiplying  this  transfer 
function  by  the  Fourier  Transform  of  a  pulse  and  obtaining  the  inverse  transform 
of  the  result. 

Figure  1  is  a  schematic  representation  of  the  measurement  scheme.  A 
Scientific-Atlanta  2020  Antenna  Analyzer  System  is  used  to  make  the  measurements. 
The  data  gathering  as  well  as  the  analysis  is  done  by  the  computer -controlled  2020 
system.  The  system  software  was  designed  in  accordance  with  Scientific-Atlanta 
Protocol  and  is  divided  into  three  sections:  1)  A  test  file  generator  section  that 
obtains  information  from  the  operator  about  the  parameters  (that  is,  position, 
carrier  frequency,  time  resolution,  etc. )  for  the  measurement  and  analysis  (Fig¬ 
ure  2  is  a  typical  listing  of  a  test  file),  2)  A  data  acquisition  section  that  acquires 
data  and  stores  it  in  a  data  file  on  disk,  and,  3)  A  data  analysis  section  that  re¬ 
trieves  the  information  from  the  data  file,  does  the  Fourier  Analysis,  and  displays 
the  results. 


(Received  for  publication  30  July  1981) 


'PULSE  RESPONSE  TEST  FILE  NR:10 

1.  POSITION 

NO  POSITION  DATA 

2.  FREQUENCY  VALUES: 

CARRIER  FREO  *  3000.00 
FREQ  INC  =  2.00 

NR  FREQ  =  512 

ACTUAL  RESOLUTION  TIME  =  .98  NANOSEC 
START  FREQ  =  2488.00 
END  FREQ  =  3510.00 

3.  RECEIVER  SETUP 

DUAL  MODE 

CRYSTAL  CURRENT  IS  PRESET 
WIDE  SEARCH 

4.  ANALYSIS 

ON  LINE  ANALYSIS  Wl  LL  BE  DONE 
FREQ  DATA  Wl  LL  NOT  BE  OUTPUT 
OUTPUT  TO  PRINTER 

SYSTEM  VALUES  WILL  NOT  BE  SUBTRACTED 
NORMALIZE  TO  CARRIER 
'•'END  OF  TEST  FILE*'* 

Figure  2.  Typical  Test  File 


2.  FOURIER  ANALYSIS 

Consider  an  ideal  pulse  p(t)  amplitude  modulating  a  carrier  at  radian  frequency 


f(t)  *  p(t)  cos  wc  t  . 

The  Fourier  transform  is 

F(w)  =  1/2  P(w  +  wc)  +  1/2  P(w  -  w )  , 

where  P(w)  is  the  Fourier  transform  of  p(t).  If  this  is  input  to  any  system  where 
the  transfer  function  is  H(w),  the  resulting  output  is 

G(w)  «  1/2  H(w)  I  P(w  +  wc)  +  P(w  -  wc)J  . 
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Hence,  in  the  time  domain 


g(t) 


l/2 


; 


j(w'  -  w  >t 

H(w'  -  wc)  P<w')  e 


dw' 


H(w'  +  wc>  P(w') 


j(w'  +  wc)t 


dw1  . 


If  we  use  the  Hermetian  properties 
H(-w)  =  H*(w) 

P(-w)  =  P*(w)  , 


we  may  write 


g(t)  « 


/ 

X 


H(w '  -  w£)  P(w') 


ejWt 


If  we  define 


A(t)  +  jB(t) 


x 

j  H(w '  -  wc)  F(w')  e'>w't  dw'  , 
-oc 


then 

g(t)  *  A(t)  cos  w  t  +  B(t)  sin  w  t 

c  c 

-  C(t)  cos  (w  t  +  ©(t)) 

c 

where. 


C(t)  *  v42(t)  +  B2(t) 


H(w'  -  wc>  P(w ') 


ejw’t 


dw' 
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and  represents  the  desired  pulse  response.  Based  upon  this  analysis  our  pro¬ 
cedure  is: 

1.  Measure  H(w) 

2.  Transform  it  to  low  frequency,  that  is,  obtain  H(w'  -  wc> 

3.  Multiply  by  the  transform  of  a  pulse,  that  is,  obtain  H(w’  -  wc>  P(w') 

4.  Take  the  inverse  Fourier  Transform  and  obtain  its  absolute  value. 


3.  MEASUREMENTS  ON  A  MICROWAVE  SYSTEM 

Figure  3  is  a  schematic  of  a  microwave  system  used  to  test  the  pulse  measure¬ 
ments  and  analysis  technique.  A  90°  hybrid  is  used  to  divide  the  signal  into:  Sj,  a 
signal  that  travels  through  a  small,  variable  and  known  time  delay,  T^  and  s2>  a 
signal  that  is  time  delayed  by  a  50  ft  cable  (time  delay  *  70  nanosec).  Also  atten¬ 
uated  by  the  cable  is  s2-  The  signals  are  recombined  in  a  magic  tee. 


Figure  3.  Schematic  of  Microwave  System  Used  to  Test  Pulse  Measurement 
System 


Figure  4  shows  a  typical  result  using  a  carrier  frequency  of  3  GHz  and  using 
512  frequency  measurements.  The  onset  of  s^  is  clearly  visible,  followed  (with  a 
70  nanosec  delay)  by  a  pulse  whose  height  is  the  vector  sum  of  Sj  and  s2>  Finally, 
the  trailing  edge  of  s2  alone  is  visible.  In  this  case  the  time  delay  T^  is  set  so 
that  the  two  signals  arrive  at  the  magic  tee  in  phase.  Thus,  the  vector  sum  equals 
the  scalar  sum.  The  signal  shown  in  Figure  4  varies  considerably  at  times  when 
smooth  "flat  top"  pulses  are  expected.  This  is  caused  by  the  measurement  system. 
Since  the  system  uses  detectors  which  are  not  ideal  and  transmission  lines  to  carry 
the  signals  to  and  from  the  detectors  the  measurement  includes  the  system  response 
that  we  call  s(w).  Thus,  the  measured  transfer  function  is 

S(w)  H(w)  . 
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Figure  4.  Pulse  Response  of  Microwave  System.  Pulse  width  =  125  nanosec, 
carrier  frequency  ■  3  GHz,  Tq  =  0 

The  system  response  S(w)  was  measured  separately,  stored  in  a  data  file,  and  used 
to  divide  out  the  system  response.  The  result  of  using  this  procedure  on  the  same 
measurement  data  that  gave  the  result  of  Figure  4  is  shown  in  Figure  5.  The  im¬ 
provement  is  clear  and  the  results  are  now  as  expected. 

Figures  5  and  7  show  results  for  values  of  T^  which  correspond  to  phase 
shifts  of  fi()°  and  120°  respectively  at  the  carrier  frequency  3  GHz.  The  effects 
of  the  vector  combination  of  the  two  pulses  is  clearly  evident.  Figure  8  shows  a 
comparison  of  the  expected  and  measured  amplitude  of  the  vector  sum  of  the  two 
pulses  as  a  function  of  the  time  delay  T^.  The  agreement  is  excellent. 


4.  ANTENNA  MEASUREMENTS 

Measurements  were  also  made  to  determine  the  pulse  response  of  two  narrow- 
band  antennas: 

1.  A  horn  and  narrow -band  filter 

2.  A  dipole  and  corner  reflector. 

Figure  9  shows  the  measured  amplitude  of  the  transfer  function  of  the  horn 
with  narrow -band  filter.  It  has  a  bandwidth  of  approximately  10  MHz  and  a  center 
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Figure  7.  Pulse  Response  of  Microwave  System.  Tq  =  l/P  nanosec 


—  COMPUTED 


0  5  1  ‘  1  -J 

0  30  60  90  120  150 

PHASE  DELAY  AT  3.0  GHz 

Figure  8.  Comparison  of  Computed  and  Measured  Vector  Sum 
of  the  Two  Pulses  Combined  by  the  Magic  Tree 
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Figure  9.  Transfer  Function  of  Horn  with 
Filter 


frequency  at  9.  38  GHz.  Figures  10  and  11  show  the  pulse  response  of  this  antenna 
for  a  carrier  frequency  of  9.  38  GHz  and  for  pulse  widths  of  125  and  250  nanosec 
respectively.  The  pulse  rise  time  is  seen  to  be  approximately  100  nanosec,  which 
is  consistent  with  the  bandwidth.  Figures  12  through  19  show  the  pulse  response 
for  various  values  of  carrier  frequency.  The  display  software  is  set  up  so  that  the 
pulse  normalization  and  pulse  width  are  variable.  In  the  above  cases  both  were 
held  constant,  the  latter  at  125  nanoaec.  Figure  20  is  a  result  at  a  carrier  fre¬ 
quency  9.  386  GHz  for  a  pulse  width  of  250  nanosec,  but  with  an  output  normaliza¬ 
tion  to  give  a  better  display  of  the  pulse  shape. 

The  measured  transfer  function  of  the  dipole  with  corner  reflector  is  shown 
in  Figure  21.  Figures  22  through  27  show  the  pulse  response  of  this  antenna  for 
various  carrier  frequencies.  The  response  for  frequencies  near  the  maximum  in 
the  frequency  response  curve  show  a  rise  time  of  about  10  nanosec  and  this  is  con¬ 
sistent  with  the  bandwidth  of  this  antenna. 

Measurements  were  also  made  on  a  wideband  antenna  (a  dipole -fed  dish  an¬ 
tenna)  but  these  results  are  not  presented  for  the  following  reasons: 

1)  At  present  the  software  only  allows  up  to  512  data  points  and  frequencies 
can  only  be  set  to  a  1  MHz  resolution.  Both  limit  the  resolution  in  time. 
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2)  The  measured  response  is  actually  the  response  of  the  transmitting  and 
receiving  antennas  combined. 

3)  Difficulties  were  encountered  in  the  technique  for  eliminating  the  measure¬ 
ment  system  response. 


TIME  {nanosec) 

Figure  10.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  *  9.  38  GHz, 
pulse  width  =  125  nanosec 


Figure  11.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  =  9.  38  GHz, 
pulse  width  =  250  nanosec 


Figure  12.  Pulse  Response  of  Horn/Filter.  Carrier  frequency*  9.  372  GHz, 
pulse  width  =  125  nanosec 
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Figure  14.  Pulse  Response  of  Horn/ Filter.  Carrier  frequency  =  9.376  GHz 
pulse  width  =125  nanosec 


Figure  15.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  =  9.378  GHz, 
pulse  width  =125  nanosec 


Figure  16.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  ■  9.  382  GHz, 
pulse  width  =125  nanosec 
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Figure  17.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  =  9.384  GHz, 
pulse  width  =12  5  nanosec 


TIME  (nanoMc) 

Figure  18.  Pulse  Response  of  Horn/ Filter.  Carrier  frequency  =  9.386  GHz, 
pulse  width  =125  nanosec 
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Figure  19.  Pulse  Response  of  Horn/Filter.  Carrier  frequency  =  9.388  GHz, 
pulse  width  =12  5  nanosec 


TIME  (nariOMc) 

Figure  20.  Pulse  Response  of  Horn/Filter  Showing  Different  Pulse  Nor¬ 
malization 
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FREQUENCY  (6  Hz) 


Figure  22.  Pulse  Response  of  Dipole  With  Corner  Reflector.  Carrier  fre 
quency  =  3.  05  GHz 


TIME  (nanosec) 

Figure  23.  Pulse  Response  of  Dipole  With  Corner  Reflector.  Carrier  fre¬ 
quency  =  3.  10  GHz 


Figure  24.  Pulse  Response  of  Dipole  With  Corner  Reflector.  Carrier  fre¬ 
quency  =  3.  15  GHz 
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Figure  25.  Pulse  Response  of  Dipole  With  Corner  Reflector.  Carrier  fre¬ 
quency  =  3.20  GHz 


TIME  (nanosec) 


Figure  26.  Pulse  Response  of  Dipole  With  Corner  Reflector.  Carrier  fre¬ 
quency  =  3.25  GHz 
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Appendix  A 

Description  of  User  Interface 

The  system  software  is  constructed  to  operate  in  accordance  with  Scientific- 
Atlanta  protocol  and  is  divided  into  three  sections. 

1.  A  test  file  generator  section  (File  No.  PUTG  in  Appendix  B)  which  obtains 
parameters  from  the  user  (antenna  position,  carrier  frequency,  time  resolution, 
etc.  )  to  be  used  for  measurement  and  analysis.  Figure  Al  shows  a  typical  test 
file  generation. 

2.  A  Data  Acquisition  Section  (File  No.  PUAQ  in  Appendix  B)  which  acquires 
data  and  stores  it  in  a  data  file  on  disk.  This  data  file  contains  system  parameters 
and  the  raw  data  —  that  is,  amplitude  and  phase  as  a  function  of  frequency. 

3.  A  Data  Analysis  Section  which  retrieves  the  information  from  the  data  file, 
does  the  Fourier  analysis,  and  displays  the  results.  (Files  No.  PUAN,  No.  PUAN2, 
No.  PUAN3,  No.  CFFT,  No.  CIFT  and  No.  PUFT  in  Appendix  B. )  Due  to  limited 
memory,  the  Analysis  Section  is  contained  in  two  segments  (second  segment  over¬ 
lays  the  first  in  memory)  with  communication  taking  place  via  system  common.  If, 
in  the  Test  File  Generation  Section,  the  user  specifies  ^n-line  analysis  then  the 
Analysis  Section  is  called  in  automatically.  However,  if  no  on-line  analysis  is 
specified,  control  passes  back  to  the  RTE  Operating  System  after  data  acquisition. 

In  this  case  analysis  can  still  be  done  by  running  PUAN  directly.  Figure  A2  shows 

a  direct  running  of  PUAN  using  Data  File  DF123  and  Figure  A3  shows  the  resulting 
output.  Caution:  In  order  to  communicate  large  blocks  of  data  between  the  seg¬ 
ments  of  the  analysis  section  it  was  necessary  to  write  over  system  common  used 


TFM  WR.PUTG 


"PULSE  RESPONSE  TEST  FILE  NR: 

1  POSITION 

WANT  TO  INPUT  POSITION  VALUES? 

VES 

TYPE  AZIMUTH  POSITION 

0 

TYPE  ELEVATION  POSITION 
30 

TYPE  AUT  PLZN  POSITION 
45 

2  FREQUENCY  VALUES 
TYPE  CARRIER  FREQ 

9500 

TYPE  REQUIRED  RESOLUTION  TIME  IN  NANOSEC 
5 

TYPE  MIN  &  MAX  FREQUENCIES 
9000,10000 

ACTUAL  RESOLUTION  TIME  =  1.95  NANOSEC 
CARRIER  FREQ  ’  9500.00 
FREQ  INC  -  100 
NR  FREQ  -  512 
START  FREQ  --  9244.00 
ENO  FREQ  »  9756.00 
3  RECEIVER  SETUP 
RECEIVER  IF  MODE 

A  *  A  ONLY,  B  “  8  ONLY,  D  =  DU/  5  =  SHARED: 
D 

RECEIVER  CRYSTAL  CURRENT  PRESET?: 

YES 


WIDE  (01  OR  NARROW  (II  SEARCH? 

0 

4  ANALYSIS 

ON  LINE  ANALYSIS? 

YES 

WANT  TO  OUTPUT  FREQ  DATA? 

NO 

RESULTS  TO  PRINTER? 

YES 

TYPE  PULSE  WIDTH  (EVEN  NO.  RES.  TIMES) 

64 

SUBTRACT  SYSTEM  VALUES? 

NO 

TYPE  1.2  OR  3  TO: 

1.  TO  NORMALIZE  TO  CARRIER  FREQ 

2.  TO  NORMALIZE  TO  MAX  AMPLITUDE 

3.  FOR  UNIT  OUTPUT  PULSE 

2 

TYPE  START  &  END  TIMES  TO  VIEW  (NO  RES.  TIMES) 

1,600 

***END  OF  TEST  FILE”' 

TFM: 


Figure  Al.  Typical  Test  File  Generation 
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PULSE  AMPLITUDE 


RU.PUAN 

TYPE  DATA  FILE  NAME 
DF123 

TYPE  LOGICAL  UNIT  NO.  (1  -  TE,  6  =  PR): 

6 

TYPE  I  TO  LIST  FREQ  DATA,  0  TO  NOT  LIST 

0 

POSITION  DATA  AVAILABLE  IN  FILE  (1  =  YES.  0  '  NO) 

0 

TYPE  1  TO  SUBTRACT  SYSTEM  VALUES,  0  TO  NOT  SUBTRACT 

0 

TYPE  1,  2  OR  3  TO- 

1.  NORMALIZE  TO  CARRIER  FREQ 

2.  NORMALIZE  TO  MAX  AMPLITUDE 

3.  UNIT  OUTPUT  PULSE 
3 

TYPE  POSITION  (NO.  RES  TIMES)  TO  NORMALIZE  TO 
256 

TYPE  START  &  END  VIEW  TIMES  (NO.  RESTIMES) 

1,500 

TYPE  PULSE  WIDTH  (NO.  RES  TIMES) 

64 

Figure  A2.  Typical  Direct  Running  of  the 
Analysis  Section 


TIME  (nanoitc) 


Figure  A3.  Output  Resulting  From  Running  the 
Analysis  Section  as  Depicted  in  A 2 
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by  Scientific-Atlanta.  At  times  this  causes  an  error  on  the  first  call  to  PUAN; 
however,  a  second  running  of  PUAN  always  gives  correct  analysis. 

The  command  files  shown  at  the  end  of  Appendix  B  are  used  to  conveniently 
integrate  the  software  into  the  2020  system.  Once  the  Fourier  programs  have 
been  compiled  and  the  relocatable  code  stored  in  files  %  PUTG,  %  PUAQ,  etc. , 
the  first  three  command  files  will  load  the  pulse  analysis  system.  The  last  two 
command  files  are  used  to  back-up  the  system  on  magnetic  tape. 


:i0 


*  4 


Appendix  B 

Software  Ustmp 


i 


#PUTG 

0001 

tttz 
0003 
0004 
0005 
0006 
0007 
0006 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
004  3 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 


T-00003  IS  ON  CR00002  USING  00021  BLKS  R-0000 

FTN4 « L 

PROGRAM  PUTG( 5  > 

C 

C  TEST  FILE  GENERATOR  FOR  ANTENNA  PULSE  RESPONSE 

C 

COMMON  I VTEK! 1 70 ) 

COMMON  IRTNt  3  )  ,  !FIG<  1  1  ,LU( 1)  ,KFNR< 1  I.IONLN!  l  > .  KDFPK! 3 > . 

1 KOF  UP ( 3  > , KOMT ( 1  )  *  I PFAL ( 1  )  , JP I A ( 1 > , IRCVRl l  >  ,KRCD! 1 >. ICMPl< 1 ). 
2 1  DUMA ( 6  > , JDISK! 1  ) . JBND< 6  )  . JBB( 6 > , JBT ( 6  > , J 1 TR ( 6  )  , JTOL ( 6  > , 
3JAZ(  1  ),UEL(  1  ) , JAUT (  1 >,JSRC( 1  >.0AUX< 1  >,JAUX2(  1  >,JSAZ( 6  >, 
4JXTAL ( 5  > . JDUM! 27  > 

C 

COMMON  I  TGN<  3  )  .  IACQ(  3  > ,  I R I D(  1  ) ,  1  U  I  DM  > ,  1 1 1  DM  > , 

1  I  2  V  0  <  1 > , I A I D<  1  ) , IE  1  D  <  1 ) , I T I D< 1  ) . 1SID< 1  ) . IUOF<  1  >. I  MODE ( 1  ). 

2  I XTAL ( 1  ), IOND( 1 > . I L I ST< 1  > , I  PLOT ( 1  > .KRMAX!  1  )  ,KUMAX<  1  ) , 

3K1MAX1  1  > , K2MAX! l  l.KAMAX! 1  J.KEMAX!  1  ),NFRO<  1  >. 

4F( 10 ) , NFR( 1  )  , I E  BM(  1  ) , I D  B  P ( l >, IDUrt< 78  ) 

C 

C 

DIMENSION  NAMSY! 3  > .KNAM! 3 ) , L NAM ( 3  > 

C 

EQUIVALENCE  ( NAMSY< 1 ) , I DUM< 30 > > 

C 

DATA  KNAM/2HPU , 2HTG , 2H  / 

DATA  LNAM/2HP  U , 2HAQ , 2H  / 

KKG-0 

C 

C  ENTRY-I FLAG  COMES  FROM  TFM 

C  IFLG-l-LIST 

C  IFLG-2-EDIT 

C  IFLG-3-WRITE 

C 

IF( IFLG  .NE .  3  )GO  TO  90 

C  SET  TEST  AND  DATA  ACQUISITION  FILES  INTO  SYSTEM  COMMON 
OO  95  1-1.3 
I TGN(  I  ) -KNAM( I  ) 

95  IACO(  M-LNAM1  I  ) 

C 

90  XGO-0 

C  FOR  FILE  CREATION  WRITE  HEADING  WITHOUT  FILE  NO. 

100  I F  <  IFLG  .EQ.  3 >WR I TE <  i . 1 0 1  ) 

C  FOR  FILE  EDIT  OR  LIST  WRITE  HEAD  I NG  WITH  FILE  NO. 

I F  <  IFLG  .EQ.  1  .OR.  IFLG  .EQ.  2  >WRITE<  1  ,  10MKFNR 

101  FORMAT! IX. "•PULSE  RESPONSE  TEST  FILE  NR:“.I2) 

C  CODE  FOR  POSITION 

KNR-1 

99  WRITE! 1 . 102 )KNR 

102  FORMAT! IX. 13, " .POSITION*  ) 

C  IF  LIST  OR  IF  ENTRY  INTO  EDIT  SIMPLY  LIST 

IF (  IFLG  .EQ.  1  .OR.  <  IFLG  .EQ.  2  .AND.  KGO  .EQ.  0MGO  TO  120 
110  WRITE! 1 . 1022) 

CALL  NOYES! JJDUM ) 

1022  FORMAT! 5X. -WANT  TO  INPUT  POSITION  VALUES?" > 

I DUM( 15 l-JODUM 
I F ( JJDUM  .EQ.  0 )GO  TO  1053 
WRITE! 1 , 105) 

105  FORMAT! 5X. “TYPE  AZIMUTH  POSITION") 

READ! 1 ,• >POSA 
KAMA X-POSA* 100.0 
WRITE! 1 , 1050) 

1050  FORMAT! 5X, "TYPE  ELEVATION  POSITION") 

READ! 1 ,*  )POSE 
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0063  KEMAX-POSEM00.0 

006*  VRITE< 1,1*51) 

0065  1051  FORMAT! 5X, "TYPE  AUT  P  LZN  POSITION" > 

0066  READ! 1 . • >POSP 

0067  KUMAX-POSPM00.0 

0068  1053  I F  <  KGO  .NE.  0)GO  TO  180 

0069  GO  TO  130 

0070  120  1FUDUMU5)  .EQ.  0)GO  TO  125 

0071  POSA-KAMAX/100.0 

0072  POSE  -K£MAX/ 100.0 

0073  POSP  -KUMAX/ 100 . 0 

0074  WRITE! l, 103 >POSA, POSE .POSP 

0075  103  FORMAT!5X, "AZIMUTH-", F8.2/5X, "ELEVATION-", F8. 2. 

0076  1/5X,"AUT  PLZN-" , F8 . 2  > 

0077  GO  TO  130 

0078  125  WRITE! 1 , 126  > 

0079  126  FORMAT! 5X, “NO  POSITION  DATA") 

0080  C  CODE  FOR  FREQUENCY  VALUES 

0081  130  KNR-KNRM 

0082  98  WRITE! 1 . 104 >KNR 

0083  104  FORMAT! IX . 13. " .FREQUENCY  VALUES:") 

0084  I F  <  I FLG  -EQ.  1  .OR.  (IFLG  .EQ.  2  .AND.  KCO  .EO.  0))GO  TO  150 

0085  140  WRITE!  1 , 106  > 

0086  106  FORMAT! 5X. "TYPE  CARRIER  FREQ") 

0087  READ! 1 , *  )FCC 

0088  WRITE! 1 , 1060) 

0089  1060  FORMAT! 5X. "TYPE  REQUIRED  RESOLUTION  TIME  IN  NANOSEC") 

0090  READ!  1 . *  )RSEC 

0091  1064  WRITS!  1 . 1061  ) 

0092  1061  FORMAT! 5X. "TYPE  MIN  &  MAX  FREQUENCIES") 

0093  READ!  1 .*  1FMIN.FMAX 

0094  NFRQ-1000. /RSEC 

0095  IDFRQ-1 

0096  I RT-NF  RQ 

0097  1062  IF! I RT  .LE.  512)GO  TO  1063 

0098  I  OF  RQ* I DF  RQ* 1 

0099  IRT-NFRO/IDFRQ 

0100  GO  TO  1062 

0101  1063  NFRQ-IRT 

0102  FDEL-FMAX-FCC 

0103  l F ! FCC-FMIN  .LT.  FDEL  >FDEL-FCC-FMIN 

0104  IFDEL-FDEL 

0105  ICO-0 

0106  ISIZ-NFRQ 

0107  1069  ISIZ-ISIZ/2 

0108  ICO-ICOM 

0109  IF! I S I Z  .NE.  1 )GO  TO  1069 

0110  NDF2-2**! ICO-1  ) 

0111  I F ! NDF  2* I  OF  RQ  .LT.  IFDEL >GO  TO  1066 

0112  WRITE! 1 . 1065 > 

0113  1065  FORMAT! IX, "RESOLUTION  TIME  TOO  SMALL  FOR  ALLOWED  FREO  RANGE") 

0114  GO  TO  140 

0115  1066  NF  RQ-2*NDF  2 

0116  1067  IF(NFRQ*IDFRQ  ,CT.  IFDEL )GO  TO  1068 

0117  I F ! NFRQ  .GT.  256  )GO  TO  1052 

0118  NFRO-2*NFRQ 

0119  GO  TO  1067 

0120  1*52  IDFRQ-IDFROM 

0121  GO  TO  1067 

0122  1068  CONTINUE 

0123  F ( 1  )  -FCC- IDFRG*NFRO/2 

*124  F!2)-FCC*IDFRQ*!NFRQ/2-l  ) 

0125  F ! 3  )- 1 OFRQ 

0126  RSEC- 1000 . /F<  3 )/NFRO 


0127 
0120 
0129 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 
01  72 
0173 
0174 
0175 
0176 
0177 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 


WRITE! 1 ,119  )RSEC 

119  FORMAT! 5X. "ACTUAL  RESOLUTION  TIME-* ,F7 .2 , ■  NANOSEC" > 

F ( 4 >-RS£C 
F(5)-FCC 
F ( 6  >*FM1N 
F ( 7 l-FMAX 
KKG-1 
GO  TO  150 

1181  IF! KGO  .NE.  0 )GO  TO  180 
GO  TO  160 

150  NFRO-! F< 2 >-F( 1 ) )/F(  3 >M .5 
WRITE! 1 , 107  )  F !  5  ) ,  F !  3  )  , NFRQ 

107  FORMAT! 5X , ‘CARRIER  FREQ«" ,F8.2/5X,'FREQ  I NC*  * ,  F8 . 2/5X , 
l  "HR  FREQ* “  ,  15  > 

FCC-F! 5  > 

WRITE! I, 1071 )F ! 1  ),F(2) 

1071  FORMAT! 5X, “START  FREQ-" ,F7 .2/5X , "END  FREO-",F7.2) 

I F  <  KKG  .EQ.  1 )GO  TO  1181 
KKF-0 

C  CODE  FOR  RECEIVER  SETUP 
160  KNR-KNRM 

WRITE! 1 , 162  )KNR 

162  FORMAT! IX, 13, -.RECEIVER  SETUP") 

IF! I F LG  .EQ.  1  .OR.  ! IFLG  .EQ.  2  .AND.  KGO  .EQ.  0>>GO  TO  165 

166  WRITE (1,1 63) 

163  FORMAT(5X, “RECEIVER  IF  MODE.-/ 

17X.-A-A  ONLY.B-B  ONL Y , D-DUAL , S- SHARE D : "  ) 

READ!  1,.  164  )IMODE 

164  FORMAT! Rl) 

IF! I MODE  .NE.  101B  .AND.  IMODE  .NE.  102B  .AND.  IMODE  .NE. 
1104B  .AND.  IMODE  .NE.  1236)166,167 

167  WRITE! 1,168) 

168  FORMAT!5X. "RECEIVER  CRYSTAL  CURRENT  PRESET?:-) 

CALL  NOYES! IXTAL  ) 

WRITE! 1 ,1680) 

1680  FORMAT! 5X, "WIDE! 0)  OR  NARROW! 1 )  SEARCH?-) 

RE AO! 1 , *  )OXTAL ( 4  ) 

I F  <  KGO  .NE.  0 )GO  TO  180 
GO  TO  190 

165  IF! IMODE  .EQ.  101 B )WR I TE < 1 , 1 65 1 ) 

I F ! IMODE  .EQ.  102B >WR I TE ! 1 . 1 652 > 

I F !  I  MODE  .EQ.  104B >WR I TE ( 1 , 1 65 3 ) 

I F ( IMODE  .EQ.  123B  )WRITE(  1  .  1654  ) 

1651  FORMAT! 8X.-A-ONLY"  ) 

1652  FORMAT! 8X. “B-ONLY"  ) 

1653  FORMAT! 8X, “OUAL  MODE") 

1654  F0RMAT(8X, “SHARED  MODE"  ) 

IF ( IXTAL  .EQ.  1 >WRITE! 1 . 1655 ) 

IF!  IXTAL  .EQ.  0 >WR I TE (  1 , 1 656  )  - 

1655  FORMAT! 8X, "CRYSTAL  CURRENT  IS  PRESET") 

1656  FORMAT! 8X, "CRYSTAL  CURRENT  NO  PRESET") 

I F ! JXTAL <  4  )  .EQ.  0 )WR I TE !  1 , 1 657 ) 

I F ( JXTAL ! 4  )  .NE.  0 )WR I TE < 1 , 1 658 ) 

1657  FORMAT18X. "WIDE  SEARCH") 

1658  FORMAT! 8X, “NARROW  SEARCH"  ) 

C  CODE  FOR  ANALYSIS  SECTION 

190  KNR*KNR*1 
WRITE! 1.191  )KNR 

191  FORMAT! IX, 13, ".ANALYSIS" ) 

IF! IFLG  .EQ.  1  .OR.  (IFLG  .EQ.  2  .AND.  KGO  .EQ.  0>>GO  TO  195 

193  WRITE! 1,194) 

194  FORMAT ! 5X , "ON  LINE  ANALYSIS?") 

CALL  NOYES! IOND) 

IANAL-IOND 


34 


0191 

0192 

0193 

0194 

0195 

1741 

0196 

0197 

0198 

0199 

1742 

0Z00 

0201 

0202 

0203 

1747 

0204 

0205 

0206 

0207 

0200 

1743 

0209 

0210 

1745 

021 1 

1  746 

0212 

198 

0213 

0214 

0215 

0216 

0217 

1744 

0218 

0219 

1952 

0220 

1951 

0221 

0222 

0223 

1953 

0224 

0225 

195 

0226 

196 

0227 

0228 

0229 

197 

0230 

1971 

0231 

0232 

0233 

1972 

0234 

1973 

0235 

0236 

0237 

1974 

0238 

0239 

0240 

1072 

0241 

1975 

0242 

0243 

0244 

0  245 

1976 

0246 

1977 

0247 

1970 

0248 

0249 

0250 

0251 

1979 

0252 

1002 

0253 

0254 

1801 

WRITE! 1 ,1741  > 

FORMAT ($X.*VANT  TO  OUTPUT  FREO  DATA?*) 

CALL  NOVES(UUDUM) 

I  0UM< 20 )"J JDUM 
WRITE! 1,1742) 

FORMAT! SR, "RESULTS  TO  PRINTER?* ) 

CALL  NOYES! JJOUM) 

IOUM! 21 1-00DUM 
WRITE! 1 , 1 747 ) 

FORMAT! 5X. "TYPE  PULSE  VI0TH<£VEN  #  RES.  TIMES)* > 

READ! 1 .*  >IPW 
IDUM! 25  >»IPW 
WRITE! 1,1743) 

FORMAT! SX," SUBTRACT  SYSTEM  VALUES?") 

CALL  NOVES(JJDUM) 

IDUM! 23  >«JJOUM 

IF !  JODUM  .EO.  0  )GO  TO  198 

WRITE! 1 ,1745) 

FORMAT!5X. "TYPE  SYSTEM  DATA  FILE  NAME  * > 

READ! I .1746  > NAMSY 
FORMAT! 3A2  ) 

WRITE!  1,1744  ) 

FORMAT! 5R. "TYPE  1,2  OR  3  TO  * " /6X . " 1 . TO  NORMALIZE  TO  CARRIER  FREO* 
1 6R . "  2 .  TO  NORMALIZE  TO  MAX  AMPL ITU0EV6X, "3. FOR  UNIT  OUTPUT  PULSE- 
READ!  1 , • >IDUM< 24  ) 

IF ! IOUM! 24  )  . NE .  3 >GO  TO  1951 
WRITE! 1,1952) 

FORMAT! 8X, "TYPE  POSITION!#  RES  TIMESJTO  NORMALIZE  TO") 

READ! I , * ) I DUM< 26) 

WRITE! 1,1953) 

FORMAT! 8R. -TYPE  START  &  END  TIMES  TO  VIEW!#  RES  TIMES)*) 

READ! 1 , * ) IDUM! 27  > , J DUM( 28  ) 

GO  TO  180 

IF! IOND  .EQ.  0 )WR I T£( 1 , 1 96 ) 

IF! IONO  .HE .  0  )WRI7E  t 1.197) 

FORMAT! 5X, -NO  ON  LINE  ANALYSIS"  i 

FORMAT! 5X, "ON  LINE  ANALYSIS  WILL  BE  DONE " > 

IF  <  IOUM!  20 )  .EQ.  1)WR I TE ! 1 . 1 97 1  > 

I F !  I OUM! 20 )  .EO.  0  >WR I TE (  1 ,  1 972  > 

F0RMAT!5X. "FREQ  DAT A  WILL  BE  OUTPUT*) 

FORMAT! 5X . "FREO  DATA  WILL  NOT  BE  OUTPUT- ) 

I F  <  I OUM! 2 1  )  .EO.  1  )  WRITE!  1,1973) 

I F ! I OUM! 21)  .EO.  0>WR5T£! 1 ,1974) 

FORMAT!  SX.'-OUTPUT  TO  PRINTER*) 

FORMAT! 5X, “OUTPUT  TO  TERMINAL") 

PUWD-IDUM! 25  J-RSEC 
WRITE! 1 .107Z iPUWO 

FORMAT!  5X.  "PULSE  WIDTH*  *  ,F8 . 2  ,  NANQSEC"  ) 

I F ( IOUM! 23  >  .EQ.  1 )WR I TE ( l , 1  975  ) 
l F <  I0UM<  23  )  .EQ.  0)WRIT£(  1 . I  976  ) 

FORMAT! 5X, "SYSTEM  VALUES  WILL  BE  SUBTRACTED*) 

FORMAT! 5X. -SYSTEM  VALUES  WILL  NOT  BE  SUBTRACTED") 

IF !  IOUM! 24  )  .EO.  1  >WR I TE l  1  ,  1 97 7  ) 

IF  <  IDUM! 24  )  .EQ,  2  )WR I TE <  1  ,  1 9 76  ) 

I F !  I OUM! Z 4  )  .EQ.  3 >WR 1 TE < 1 , 1 979  ) 

FORMAT (5X.-HORMAU2E  TO  CARRIER") 

FORMAT($X, "NORMALIZE  TO  MAX  AMPLITUDE") 

FORMAT! SX. "UNIT  OUTPUT  PULSE") 

I F ( I OUM! 2 4  )  . N£ .  3  )GO  TO  1801 
TN0R*RS£C*1DUM(26) 

WRITE! 1 , 1802 )TNOR 

FORMAT! 8X, "NORMAL IZE  TO  OUTPUT  AT*.F8.2."  NANOSEC") 

ST IM»RS£C* 1 DUM( 27  ) 

ET1M-RSEC*1DUM( 28  ) 
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1255 

WRITE! l . 1003 )ST1M,ETIM 

#256 

1003 

1  FORMAT(  8X,  "VIEW  START  TIME-*  ,F8. 2  .  'NANOSECV 

#257 

18X.-VIEW  END  TIME-* ,F8.2.“  NANOSEC") 

0258 

GO  TO  180 

0259 

180 

WRITE( 1.108) 

0260 

108 

FORMAT! 1 X. "•••END  OF  TEST  FILE****) 

0261 

GO  TO<  250. 1 70 , 250  )  I F LG 

0262 

170 

WRITE! 1 . 109 ) 

0263 

109 

FORMAT! IX. -TVPE  LINE  NO.  TO  EDIT  (-1  TO  END)' 

0264 

READ! 1 ,• ) ICC 

0265 

KGO-1 

0266 

ICC- ICC* 1 

0267 

GO  TO!  250. 110.140, 166. 193. 250)  ICC 

0266 

C 

RETURN  TO  TFM 

0269 

250 

CALL  EXEC! 8 , IRTN  ) 

02  70 

CALL  MAIN 

0271 

END 

0272 

ENDS 

:i6 


>  -j 


##  63  C 

#064  C  BEGIN  FREQ 

##65  C 

##66  MDFMt 33  >-2HBF 

##67  MDFMt  35 ) *32 

##68  C 

##69  C  E NO  FREQ 

887 B  C 

#071  MDFMt  39  )*2MEF 

##72  MOFM(41)-32 

##73  C 

##74  C  STEP  FREQ 

##75  C 

##76  MDFM( 45 ) ■2HSF 

##77  MDFM(47)-32 

0078  C 

##79  C  DATA  POINT  MARKER 

aaaa  r 

#081  MDFM(49 J-NFRQ 

##82  MDFMt 51  1-2H0P 

0083  MDFMt 52)-# 

#08  4  C 

0085  C  FREQ 

0086  C 

0087  MDFM<54>--2 

0088  MDFMt 56 >-2HFR 

0089  MDFMt  58  >■! 6 

009#  C 

0091  C  AMPLITUDE 

#092  C 

#093  MDFMt  6#)«-l 

0094  MDFMt 62 )-ZHAA 

0095  MDFMt 631-3 

0096  MDFMt  64  ) ■- 1 4 

009  7  C 

0098  C  PHASE 

0099  C 

0100  MDFMt  68 )•- J 

0101  MDFMt 701-ZHPA 

0102  MDFMt  711-3 

#103  MDFMt  72  )*-!3 

#104  C 

0106  C  SET  DISK 

0106  C 

#107  JDISK-1 

0108  C 

#109  C  TEMP  SET  FOR  NO  BMAX 

#1  18  C 

#111  MDFMt 14  )■# 

0112  C 

0113  C  SET  TF  FILE 

0114  C 

0115  ITGNt  D-2HPU 

#116  ITGN(2)-2HTG 

#117  ITGNt  3  )»2H 

#118  ITGN(4)-2HPU 

#119  ITGN(S)-2HAQ 

#12#  ITGNt  6  1-2H 

#121  C 

0122  C  OPEN  FILE  6  TRANSFER  FORMATS 

#123  C 

#124  1#  CALL  DFPKt 1 , 1 , 1 ,LERR  ) 

#125  C 

#126  C  TRANSFER  TF 


#PUAQ  T-«f««3  IS  ON  CR00002  USING  09921  BLKS  H-fMf 


fill 

1112 
1113 
0004 
0005 
0006 
000 7 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 


FTN4.L 

PROGRAM  PUAQ(S) 

COMMON  I VTE  K( 170) 

COMMON  I RTN( 3  > ,  I FLG< 1 > ,LU( 1 ) , KFNR< 1 > . XONLNC 1 ),KDFPK<  3>. 

1 KDF  UP ( 3  > ♦ KOMT ( l  >,IPFAL< 1 J.JPIAI 1 ) , IRCVR( 1 ) ,  KRCDI 1 ) . 1CMPLI 1 1 , 
21 0UMA( 6 ) ,  JDISKI 1  ) , JBNDI 6  )  ,0BB< 6  )  .OBTI 6  )  ,JITR<6  >  ,JTOL<  6  > . 

3  JAZ ( 1 >.JEL(  1 ) , JAUT ( 1 >,J$RC( 1 > ,JAUX< 1 >.0AUX2< 1  >,JSAZ<  6), 

4 JXTAL ( 5  > .  JDUMI 27  > 

COMMON  ITGN( 3  )« IACQ< 3), IRID( 1 ) , IU ID( 1 ) . I  1 ID<  1  ) , 

1I2I0U  >.  IAID<  1  >,  IE  10(1  > , I TIDI  1  ) ,  I S 1 0  f  1  >.IU0F<1  ).IMOOE(1  ). 

2 1 XTAL ( 11. IOND< 1 >, ILISTI 1  )  ,1P10T( 1 > , KRMAX( 1 >,KUMAX< 1 ) , 

3K1MAXI 1 >.K2MAX( 1 ) , KAMAX1 1 ),KEMAX( 1 >.NFRO( 1 ). 

4F< 10) , NFR( 1 ) , I EBM< 1 ).I08P<  2 ).IDUM(78) 

COMMON  IEEE( 1 28  ) , MIAMI  8  )  . MDFMl 120) 

DIMENSION  KDNAM( 3  )  • KTFM< 3 ) , KANAMi 3 ) 

DIMENSION  FF(3).ICCC 128), I OAK  8  )  ,KKNAM( 5 ) , POS<  3 ) 

EQUIVALENCE  ( I A , I  OAT  <  3  )  ) » ( I MA . I DAT ( 4  )  ) , ( ISA, I DATC  5  > ) 
EQUIVALENCE  (  IP  .  I  OAK  6  )  > .  <  I  PR ,  I  DAK  7  >  ) ,  (  I L  S .  IDAT<  8  )  > 

DATA  KKNAM/2HCO , 2HMM , 2HT  ,2HTF,0/ 

DATA  KDNAM/2HDF , 2HM  , 2H  / 

DATA  KTFM/2HTF , 2HM  ,2H  / 

DATA  KANAM/2HPU , 2HAN , 2H  / 

C 

C  INITIALIZE  FORMATS 
C 

DO  1  1-1.8 

1  MIAMI  I )-0 

DO  2  1-1,120 

2  MDFMl I )-0 
DO  3  1-1.74 

3  MDFMl  I  )•»  1 
C 

C  TEST  FILE 
C 

MDFMl  2)  — 128 
MDFM( 4  )-2HTF 
MDFM( 6 >-400008-15 

C 

C  COMMENTS 

C 

MDFMI8)— 128 
M0FM< I01-2HCO 
MDFM( 12  )-16 

I F C KOMT  .EQ.  0)MDFM<8l-0 
C 

C  BEAM  MAX 

C 

MDFMl  14)  — 128 
MDFMl 161-2H8M 
MDFM< 18 )-40000B-15 
C 

C  POSITION  MARKER 
C 

MDFM<  22  )-2HPM 
MDFMl 23 )-0 
C 

C  POSITION 
C 

MDFMl  25  >  —  3 
MDFMl 27  )-2HPO 
MDFMl 29)-32 

IFIIDUMI1S)  .EQ.  0  )MDFM<  25 )"0 
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0127 

C 

0128 

call  DFP  K  ( 2, ITGN.128.LERR) 

0129 

C 

0130 

c 

TRANSFER  COMMENTS 

0131 

c 

0132 

IFtKOMT  .EQ.  0>GO  TO  45 

0133 

CALL  DFPKt 3.KKNAM, 128, LERR ) 

0134 

45 

IFUOUMU5)  . EQ .  0>GO  TO  46 

0135 

C 

0136 

C 

TRANSFER  POSITION 

0137 

C 

0138 

POS( l ) -KAMAX / I 00 . 0 

0139 

POS( 2 l-KEMAX/ 100.0 

0140 

POS( 3  )-KUMAX/ 100.0 

0141 

DO  452  1-1.3 

0142 

ANG*PO$< l > 

0143 

IAX*I 

0144 

CALL  PSNt 1AX.ANG, ISTAT) 

0145 

I F <  ISTAT  .NE .  0>VRITE< 1,451 ) 

0146 

CALL  SVN0( ANGM. IAST) 

0147 

CALL  SYN1 

0148 

I F ( IAST  .NE.  0)WRITE( 1 ,451 > 

0149 

POS< I 1-ANGM 

0150 

452 

CONTINUE 

0151 

451 

FORMATt IX, -POSITION  ERROR") 

0152 

CALL  DFPKt 2.POS.6.1ERR) 

0153 

C 

0154 

C 

TRANSFER  BEG. END, STEP  FREQ 

0155 

C 

0156 

46 

DO  47  1*1,3 

0157 

47 

F F (  I  ) *F <  I  ) 

0158 

CALL  DFPKt 2.FF ,6,LERR ) 

0159 

C 

SET  VARIABLES  FOR  LOWEST  ATTENUATION 

0160 

I  FA-0 

0161 

IFB-0 

0162 

C 

0163 

c 

SET  UP  &  TRANSFER  F RE Q . AMP , P HASE  DATA 

0164 

c 

0165 

c 

SET  UP  RECEIVER 

0166 

c 

CHOOSE  REMOTE  TUNING 

0167 

ILR-I 

0168 

c 

ENABLE  APC  SEARCH 

0169 

IAPC-1 

0170 

c 

USE  WIDE  OR  NARROW  SEARCH  DEPENDING  ON  JXTAI41 

0171 

c 

DISABLE  REFERENCE  CHECK 

01  72 

IREF-0 

01  73 

c 

CHOOSE  MODE 

0174 

IFF--1 

0175 

I  F  (  I  MODE  .EQ.  101BHFF.J 

0176 

IFtIMODE  .EQ.  102BHFF-2 

0177 

I  F  t  I  MODE  .EQ.  104BUFF-3 

0178 

IFtIMODE  .EQ.  123B  > IFF-0 

0179 

CALL  RMOOEt ILR . IAPC .OXTALl A > , IXTAL . 1REF , IFF : 

0180 

c 

SET  RECEIVER  TO  AUTO 

0181 

CALL  RMANt 1 ) 

0182 

FZ-FF( 1 > 

0183 

CALL  FR£Q<0.FZ.FCNTR, ISTAT) 

0184 

c 

SET  OATA  TO  1 

0185 

DO  11  1-1,8 

0186 

11 

1  OAT  t  1  )-l 

0187 

c 

INITIALIZE  RATIOMETER 

0188 

CALL  RAT0< IA, IMA. ISA, IB.IMB.ISB.ISTA) 

0189 

c 

INITIALIZE  PHASE  DISPLAY 

0190 

CALL  PHA0( IP. I PR . IIS, ISTP ) 
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*191 

*192 

*193 

*194 

*195 

*196 

*197 

*198 

*199 

*20* 

*2*1 

*2*2 

*2*3 

*2*4 

*2*5 

*2*6 

*2*7 

*2*8 

*2*9 

*210 

*211 

*212 

*213 

*214 

*215 

*216 

*217 

0218 

*219 

*22* 

*221 

*222 

*223 

0224 

*225 

*226 

*227 

*228 

*229 

*23* 

*231 

*232 

*233 

*234 

*235 

0236 

*237 

*238 

*239 

*24* 

*241 

0242 

*243 

*244 

*246 

*246 

*247 

*248 

*249 

*25* 

*251 

*252 

*253 

*254 


C  DO  MEASUREMENTS  FOR  EACH  FREQUENCY 
PCTOL "F ( 4 ) 

DO  12  1-1  ,  NFRQ 
FZ-FF<  1  )♦<  1-1  )*FF( 3) 

C  SET  FREQUENCY  AND  MEASURE  IT 
1TRY-* 

5*  CALL  F REQ( *, FZ, FCNTR , ISTAT ) 
IFZ-FZ 

IFCN*FCNTR*.5 

IF! I FZ  .EQ.  IFCN1G0  TO  53 

ITRY-ITRVM 

I F ( I TRY  .LT.  l**)GO  TO  5* 

53  CONTINUE 
C  RECORD  MEASURED  VALUE 
I  DAT (  1  )  -FCNTR 

IDAT!2>-32767*!FCNTR-IDAT! 1>> 

C  SET  ATTENUATION 

CALL  RATN! IFA, IFB ) 

CALL  DLAY! 1**> 

C  CHECK  RECEIVER  STATUS 

CALL  RSTAT! ILOK, IORNG) 

IJU-ILOKM 
GO  TO< 1*0 , 288 ) 1 0 U 
10*  WRITE! 1,15) 

15  FORMAT! IX, “NO  PHASE  LOCK*  ) 

2*0  CONTINUE 
C  TRIGGER  RATIOMETER 
29*  CALL  RATI 
C  GET  DATA 
I  TRY"* 


30*  CALL  RAT2 

IF!  ISTA  .EO.  * )GO  TO  55* 

I F <  ITRY  .EQ.  1*1  )GO  TO  54* 

CALL  DLAY ( 10) 

ITRY" ITRV+1 
GO  TO  3*0 

54*  WRITE! 1 ,50*0) 

500*  FORMAT! IX. “RATIO  TIMEOUT"  ) 

GO  TO  1*0* 

C  DB  VALUES  MUST  BE  NEGATIVE 
55*  IF! IB  .GE.  *>IFB-IFB*t 
IF! I A  .GE.  *  > I  FA" I  FA* 1 

I F !  I A  .LT.  *  .AND.  IB  .LT.  *)CO  TO  58* 
IF!  IFA  .EQ.  7  .OR.  IFB  .EQ.  7 >GO  TO  59* 
CALL  RATN! IFA, IFB  ) 

CALL  OLAY! 1**) 

GO  TO  290 

59*  WRITE! 1 ,55**) 

550*  FORMAT! IX, "AMPLITUDE  TOO  HIGH") 

GO  TO  1*0* 

58*  CONTINUE 
1ATE-IA 
IBTE-IB 
IA-IA45**IFA 
IB-IB^5**IFB 
IF1IATE  .LT.  -3**>IFA«* 

IF ( I BTE  .LT.  -3**)!FB«* 

C  TRIGGER  PHASE  DISPLAY 
CALL  PHA1 
C  GET  DATA 


ITRY-* 

C  CHECK  STATUS  OF  PHASE  METER 
45*  CALL  PHA2 
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0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 


I F (  ISTP  .£Q.  0X50  TO  480 
I F ( ITRY  .GE .  101 >C0  TO  470 
CALL  DLAV< 10) 

ITRV" I  TRY* 1 
GO  TO  450 

470  WRITE! 1 .4500) 

4500  FORMAT! IX. -PHASE  TIMEOUT-) 

GO  TO  1000 
480  CONTINUE 
C  TRANSFER  DATA 

CALL  DFPK! 2 , I  DAT , 8 , LERR  > 

12  CONTINUE 
C  CLOSE  DATA  FILE 

CALL  DFPK!  8 , 1  , 1 , LERR  ) 

ICMPL-l 

I F ( I ONO  .NE.  0)400.500 

1000  WRITE! 1 . 1001 ) 

1001  FORMAT! IX. "DATA  FILE  ABORTED- > 

CALL  DFPK! 7, iDUMY, IDUMV.LERR) 

WRITE! 1 , 1002 ) 

1002  FORMAT!  IX. "TRY  AGAIN?- ) 

CALL  NOYES! ITRY  > 

IF! ITRY  .NE .  0  )G0  TO  10 
500  CALL  EXEC! 8 , KTF M  ) 

CALL  MAIN 

400  KDFUP!  1  >"KDFPK! 1  ) 

KDFUP ( 2  >*KDFPK( 2  > 

KDFUP! 3  J-KDFPK! 3  > 

M»0 

1005  I CODE ■ 1 00000B ♦ 9 

CALL  EXEC! ICODE .KANAM) 

M"  1 

I F  <  M  .EQ.  1 >G0  TO  405 
GO  TO  500 
405  WR I TE ! 1 . 401 ) 

401  FORMAT! IX, -ERROR  IN  SCHEDULING  ANALYSIS  PROGRAM-. 
1/2X, "TRY  AGAIN?-  ) 

CALL  NOYES! ITRY) 

IF!  ITRV  .NE .  0 )G0  TO  1005 

GO  TO  500 

END 

ENDS 


41 


#PUAN  T *00003  IS  ON  CR0000?  USING  0002 1  BLKS  R-flff 


mi 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  l 
0012 
0013 
0014 
00!  5 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
002  7 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
004  4 
004  5 
0046 
004  7 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 


FTNi.L 

PROGRAM  PUANO) 

COMMON  I VTEK! 1 70 > 

COMMON  1 RTN( 3  > , 1 FLC ( 1  ) . LU! 1  ) . KF NR( 1  ). IONLN( 1 >,K0FPK!  3>. 
1KDFUP(3), KOMT  < I  >,IPFAL<  I J.JPIA! 1 ),IRCVR! 1 >,KRCD( 1 >,ICMPL< 1 >. 
2I0UMA(6  ). JDISKt 1  > . JBND! 6 > . JB8( 6 ) ,0BT! 6  ) , JITR! 6 ) ,JTOl! 6  ) . 

3JAZ!  1  > .  OEL  <  1  )  ,  JAUT!  1  ) .  JSRCM  )  ,  JAUX!  1  )  ,  JAUX2!  1  > ,  JSAZ!  6  ) , 

4JXTAL ( 5  > , JDUMI 27  ) 

COMMON  ITGN( 3 ) , 1ACQ<  3  > . 1 R I D( 1 > • I U1 DC  1  ) . 1110(1). 

1 12I0< It, IAIO< 1 J.IEID! 1 >. ITIO! 1 >, I  SID!  I >. IUOF! 1  ). IMODE! 1 >. 

2 I XTAL <  1  )  . I ONO (  1  >. 11 IST< 1  >. 1PL0T( 1  )  .KRMAX!  1  > .KUMAX! 1 >, 

3K 1  MAX (  1  >,K2MAX( 1  ),KAMAX( 1  ),KEMAX( 1  > , NFRQ<  1  ). 

4F( 10 ) , NFR! 1  )  , 1 EBM(  1 > , I D8P ( 1 > , 1 0UM< 78 ) 

COMMON  ! FEE ( 128 ) .MIAMI  8 > ,MDFM( 120) 

0 IME  NS  ION  DB(  51 2  >.  DEG!  51 2  > 

DIMENSION  ICC! 128).FFA! 3  ) ,  I  DAT ( 8  > , POS ( 3  ) 

DIMENSION  KANAM! 3  ) 

COMPLEX  X<  51 2  > 

EQUIVALENCE  ( !CC< 1  ) , MIAMI  1  ) ) 

DATA  KANAM/2HPU, 2 HAN  .  2H 1  / 

PI-3.1415926 

JSKIP-0 

I F  <  I OND  .NE .  0)4,3 
C  OBTAIN  DATA  FILE  NAME 
3  WRITE! 4,60) 

60  FORMAT! IX. "TYPE  OATA  FILE  NAME') 

READ! I . 70JKDFUP 

70  FORMAT! 3A2) 

C  OUTPUT  TO  TERMINAL  OR  PRINTER? 

WRITE! 1,5) 

5  FORMAT! IX, 'TYPE  LOGICAL  UNIT  NO. ( 1-TE ,6-PR  ) : *  ) 

READ!  J  ,*UUA 
IDUM{ 22 )-LUA 
WRITE! I .71  ) 

71  FORMAT! IX, "TYPE  1  TO  LIST  FREQ  DATA , 0  TO  NOT  LIST") 

READ! 1 ,*)JO20 

WRITE!  1  ,72  ) 

72  FORMAT! IX, -POSITION  DATA  AVAILABLE  IN  FILE! 1-YES  ,0-NO )*  ) 

READ! 1 ,•>0015 

WRITE! 1 .73  ) 

73  FORMAT! IX, "TYPE  l  TO  SUBTRACT  SYSTEM  VALUES.0  TO  NOT  SUBTRACT*) 
READ! I .*  )0D23 

IF ! JD23  .EQ.  0 )G0  TO  75 
WRITE! 1 .76) 

76  FORMAT! IX. "TYPE  SYSTEM  OATA  FILE  NAME") 

READ! 1 , 77  )IDUM! 30) , ZDUM! 31 >, I  DOM! 32  > 

77  FORMAT! 3A2 ) 

75  WR 1 T E  < 1 , 7 4  ) 

74  FORMAT! IX. "TYPE  1.2  OR  3  TO- " / 3X , " l . NORMAL  I ZE  TO  CARRIER  FREQ"/ 
13X, "2. NORMALIZE  TO  MAX  AMPL I TUDE " /3X , " 3 . UN  I T  OUTPUT  PULSE") 

READ! 1 , • >0024 
I DUM<  24 ) ■ J024 
IF  <  JD24  .NE.  3  )C0  TO  78 
WRITE! 1,79) 

79  FORMAT! IX, "TYPE  POSITION!#  RES  TIMES  )T0  NORMALIZE  TO") 

READ! 1 («>1DUM!26> 

78  WR I TE ( 1 , 79 1 ) 

791  FORMAT! IX, -TYPE  START  4  END  VIEW  TIMES!#  RES  TIMES)") 

READ!  1 ,*  )  IDUM! 27  ), IDUMl 28  > 

WRITE! 1,792) 

792  FORMAT! IX . "TYPE  PULSE  WIDTH!#  RES  TIMES)") 

REAO! 1 . • >IDUM< 25 ) 


42 


t 

k 


1163 

GO  TO  6 

1164 

4 

IOND-0 

0065 

J  D 1  5  •  1 DUM!  15  ) 

0066 

JO20-IDUH( 20) 

006  7 

J023-1DUM! 23  ) 

0068 

J D24 ■ I DUM! 24  > 

0069 

LUA-1 

0070 

I F  < I DUH<  21)  .EQ.  1  )LUA»6 

0071 

1 0UM( 22 )-LUA 

0072 

6 

CONTINUE 

0073 

C 

OPEN  FILE  &  GET  FORMATS 

0074 

CALL  DF UP ( I , LERR  > 

0075 

KERR-1 

0076 

1F( LERR  .HE .  0>GO  TO  101 

0077 

C 

READ  4  OUTPUT  TEST  FILE 

0078 

CALL  DFUP! 4, LERR, ICC. 128) 

0079 

WR  I TE ( LUA , 10  )  ICC 

0080 

10 

FORMAT! 1 X . 32A2 > 

0081 

KERR-2 

0082 

1 F ( L E RR  .NE.  0)GO  TO  101 

0083 

c 

READ  &  OUTPUT  COMMENTS 

0084 

IF<KOMT  .EQ.  0 )GO  TO  35 

0085 

CALL  DFUP! 4.LERR, ICC. 128) 

0086 

WRITE! LUA, 10)ICC 

0087 

KERR-3 

0088 

IF ( LERR  .NE.  0)GO  TO  101 

0089 

35 

IFIJD15  .EQ.  0  )GO  TO  25 

0090 

C 

READ  A  OUJPUT  POSITION 

0091 

CALL  DFUP! 4.LERR.POS.6) 

0092 

WRITE!LUA.20)POS 

0093 

20 

FORMAT! 10X.F8.2  > 

0094 

KERR-4 

0095 

IF(LERR  .NE.  0 )GO  TO  101 

0096 

C 

READ  START , END  A  FREQ  INC  A  OUTPUT 

0097 

25 

CALL  DFUP( 4.LERR.FFA.6) 

0098 

WR I TE ! L  UA , 20 ) F  F A 

0099 

KERR-5 

0100 

IF! LERR  .NE.  0)GO  TO  101 

0101 

NFRQ-!FFA<2  )-FFA( 1 ) >/FFA(3)M .5 

0102 

N-NFRQ 

0103 

NCO-0 

0104 

100 

M-N 

0105 

N-N/2 

0106 

NCO-NCO+1 

0107 

I F ( N*2  .NE.  M )GO  TO  105 

0108 

I F ( N  .NE .  1  )CO  TO  100 

0109 

IF! JSKIP  .EO.  1 )GO  TO  202 

01  10 

C 

FOR  EACH  FREQ  READ  A  OUTPUT: 

01  1  1 

C 

! 1  )FREO  >  IDAT! 1  ) , 1  DAT ( 2  ) 

0112 

C 

! 2  )AMP  >  IDAT! 3-5) 

0113 

C 

( 3  )PHASE  >  I  DAT ( 6-8  ) 

0114 

I F ! JD20  .EO.  1 )WRITE<LUA.45) 

0115 

45 

FORMAT! 10X, -FREQ" , 1 4X , "AMP ! OB  )  " , 1 

01  16 

202 

D8MAX-- 100 

0117 

DO  12  I-l.NFRO 

0118 

CALL  DFUP(4,LERR, I  DAT, 8  ) 

0119 

AMULT-1 .0 

0120 

I F  <  1  DAT  <  4  )  .  EQ .  DAMULT-.l 

0121 

I F  < 1  DAT  <  4  )  .EQ.  0  JAMULT- .01 

0122 

IF! AMULT  .EO.  l.0)GO  TO  101 

0123 

D8SVS- 1  DAT ! 3 )*AMULT 

0124 

DCSVS-0. 1  *  I  DAT ! 6  ) 

0125 

IF! JSKIP  .EO.  1  )GO  TO  61 

0126 

D8( I  )-DBS VS 
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*127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0179 

0179 

0180 

0181 


DEG! I  )-DGSYS 
GO  TO  62 

61  DB<  I  >-DB<  Il-DBSVS 
DEG <  I >-DEG( I  ) -DGS VS 

62  1F(DB<1)  -GT.  DBMAX  ) D0MAX -DB< I  ) 

FCC-  I  DAT  <  1 >*IDAT<  2J/32767.0 

I F  <  J  D  20  .EQ.  1  )WRITE(lUA,40)FCC,DB< I  ) , DEG ( !  ) 

40  FORMAT! 3< 10X.F8.2 > > 

12  CONTINUE 

KERR-6 

1FUERR  .NE.  0  )GO  TO  101 
CALL  DFUP!  6 , L ERR  ) 

If ( JD23  .NE .  1  >GO  TO  230 
IF( JSKIP  .EQ.  1 >GO  TO  230 
DO  201  1-1.3 

201  KDFUP< I >»IDUH( 29*1  ) 

JSKIP- l 
GO  TO  6 
230  CONTINUE 

NFRQ2-NFRQ/2 

IF  <  JD24  .EQ.  1 )DBMAX*0B( NFRQ2+1 ) 

DO  14  I-l.NFRQ 
0B<  I  )  -DB! I  ) -DBMAX 

14  CONTINUE 
DEG1-DEG!  NFRQ2M  ) 

DO  15  1-1 , NFRQ 
DEG  <  I  >-DEG< I > -DEG  I 

21  I F ! DEG( I  )  .LE.  -180.  >DEG<  I  >-DEG!  I  >060. 

IF ( DEGll  )  .GT.  180.  >DEG( I )-DEG( I >-360. 

IF! DEG< I  )  .LE.  -180.  .OR.  DEG( I  )  .GT.  180.  )GO  TO  21 

15  CONTINUE 
IWP- I DUM( 25 ) 

CALL  PUFT! X . NCO , IWP ) 

DO  19  1-1 .NFRQ 
RAD- DEG ( I >*PI/180. 

AMP-10. 0**<  DB(  M/20.0) 

IF!  I.LE.NFRQ2)X(I ♦NFRQ2 )-AMP*CMPLX( COS( RAD  ) , -SIN! RAO ) >*X!  I ♦ NFRQ 2 ) 
IF < I .GT . NF RQ2  >X! I-NFRQ2 )*AMP*CMPLX( COS( RAD  > , -SIN! RAD ) )*X( I -NFR02  > 
19  CONTINUE 

GO  TO  107 
105  WRITE! 1,102) 

102  FORMAT! IX, "NO.  OF  FREQ  NOT  A  POWER  OF  TWO- > 

GO  TO  99 

101  WRITE! 1 ,30>KERR,LERR 

30  FORMAT! IX, "ERROR  IN  DFUP  #  -tI2,"  LERR-",I5) 

99  CALL  DFUP ! 6 , LERR  ) 

108  CALL  EXEC! 6  ) 

107  IDUM! 1  )-ICO 

IOUM! 2  >«NCO 
IDUM! 3  >-LUA 
F ! 3  )-FFA( 3  > 

CALL  PUAN2! X ) 

GO  TO  108 

END 

ENDS 
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0001  FTN4.L 

0002  SUBROUTINE  PUAN2IX) 

0003  COMMON  IVTEKU70) 

0004  COMMON  I RTN ( 3), I  FLG<  1  > » L U<  1 > , KFNR< I  ). I ONLN(  1  ) ,  KDF  PK( 3  > . 

0005  l  KDF  UP  (  3  > .  KOMT  (  1  ).  IPFAL1  1  )  .  JPIAI  1) ,  IRCVR<  1  >.KRCD(  1  ).  ICMPLC  1  ), 

0006  2I0UMA( 6  )  , JDISKI  1  > . J BND < 6  )  . J 68 1 6  ) , J B T ( 6 > , J I TR ( 6  ) , J TOL ' 6  ) , 

0007  3  JAZ<  U.OEU  1  l.JAUTt  1  >.OSRC<  1  >.JAUX<  1  J.0AUX21  1  J.0SAZ16  >. 

0008  4JXTAK  5  > ,  JDUMt  27 ) 

0009  C 

0010  COMMON  ITGN< 3  >. IACQI 3 >. I R T D ( 1 > , IUID< 1  ) , 1 1 ID( 1 >, 

0011  1 1 2  I  D<  1  >  .  I A  1 0 (  1  >  .  I  E  I  D(  1  > . I T 1 D  <  1  )  , I S I D(  1  ) . I UOF (  1  ) . JMODE(  1  ) . 

0012  21XTAI (  1  ) ,  1  OND(  D.ILlSTU  >,1PL0T<  1  ),KRMAX<  1  J.KUMAXt  1  ), 

0013  3KIMAXI  I  ) , K2MAX (  1  ),KAMAX( 1  >,KEMAX< 1  > . NF  RQ(  1  ). 

0014  4F(  10  )  , NF  R(  1  ) , IEBM<  1  ). IDBP(  1  )  , IDUM< 78  ) 

0015  COMMON  IE£E<  120 ) .MIAMI  8  )  ,MDFM<  120) 

0016  DIMENSION  I  DAT (512), KANAM ( 3  ) 

0017  COMPLEX  X ( 5 1 2  ) 

0018  EQUIVALENCE  ( RSEC . JDUM( 25 > > 

0019  EQUIVALENCE  ( ITGNI  1 > . I  DAK  1) ) 

0020  DATA  KANAM/ 2 HP  U , 2HAN . 2H3  / 

0021  RSEC-1000. /( NFRQ-1  )/F< 3  ) 

0022  JDUMf 27  )-I0UM( 28  )-IDUM( 27  ) 

0023  ICO-IDUM(I) 

0024  NCO* I DUM<  2  ) 

.*025  LUA-IDUM(22> 

0026  JD24*IDUM( 24  ) 

0027  NFRQ2-NFRQ/2 

0020  CALL  CIFT(X.NCO) 

0029  XR  =  REAL( X( NFRQ2  )  ) 

0030  X  I *A IMAG ( X<  NF  RQ2  )  > 

0031  XNOR-SQRT< XR*XR*XI*XI> 

0032  IF1JD24  .NE.  3)XNOR»l. 

0033  I S* I DUM( 27 ) 

0034  I E  *  I DUM( 28  ) 

0035  DO  1  I  *  I S  ,  IE 

0036  XR«REAUX<  1  >> 

0037  XI-AIMAG1 X( I ) ) 

0038  XR»$QRT<XR*XR*XI*XI  l/XNOR 

0039  IOAT(  I-ISM  )-300.*XR 

0040  1  CONTINUE 

0041  I  CODE ■ 1000008*9 

0042  CALL  EXEC( ICODE .KANAM) 

0043  RETURN 

0044  END 

0045  ENDS 


#PUAN3  T>f|f|3  IS  OR  CR00002  USING  00008  BLKS  l-flff 


fill 

FTN4 

.1 

10*2 

PROGRAM  PUAN3I3) 

0003 

COMMON  IVTEKII70) 

0004 

COMMON  1RTN<3>,  IFLC!  t> ,  LUI  1  ) ,  KFNR(  1  KIONLN!  1  ) .  KDFPK!  3  ) , 

0005 

lKDFUP!  3  )  ,  KOMT!  1  ) .  I  PFALI  1  ) ,  JP I  A!  1  )  ,  IRCVRI  1  I.KRCDI  1  ).  ICMPU  1  > 

0006 

2IDUMA!6).J0ISKf 1  ),JBND(6>.JBB(6),JBT(6),JITR(6),JTOL(6), 

0007 

3JAZ<  1  >,JEL<  1  l.JAUT!  1  >,JSRC<  1  ),0AUX<  1  I.JAUX2!  I  ),JSAZ(6), 

0000 

40XTAL! 5 ) , JDUM1  27  > 

0009 

C 

0010 

COMMON  1TGN<3).IAC0(3),IRID(1>,IUID(1).MIDU). 

0011 

1 12I0(  1  l.IAIO!  1).  IE10<  l  ) .  ITI0(  l )  .IS  ID!  t  )  ,IUOF(  l ) ,  I MODE <  1 ). 

0012 

2IXTAL! 1  ),  IONO<  l  ),  ILIST!  I  ),  I  PLOT!  1  I.KRMAXI  1  I.KUMAXi  1  ), 

0013 

3X1 MAXI 1 I.K2MAXI 1 ) ,XAMAX( I ).KEMAX< 1 >,NFRQ( 1 ), 

001  4 

4F! 10I.NFR! 1 ), IEBM!  1). 1DBP(  1  ), IDUM( 78) 

0015 

COMMON  IEEE( 12B ) ,MIAM<  8 ) ,MOFM( 120) 

0016 

DIMENSION  I VL ( 1*) , IXA<  7 ) , IVA! 15 ) , II ( 2 > . I DAT! S 1 2 ) 

0017 

EQUIVALENCE  < I TCNI l ) , IDAT( 1) ) 

0018 

EQUIVALENCE  < RSEC . ODUM! 25  )  > 

0019 

DATA  IVA/2H  P,2H  U,ZH  L , 2H  S,ZH  E , 2H  ,2H  A,2M  M.2H  P.2M  L. 

0020 

12H  I , 2b  T,2H  u , 2H  D,2H  E/ 

0021 

DATA  IXA/ZHTI^HME.ZHIN.EHAN.ZHOS.ZHEC^H)  7 

0022 

DATA  IVL/2M0. ,2H0  .2H0..2H5  .ZH1..2H0  .2H1..2H5  .ZH2..2M0  / 

0023 

ID-J0UM(27) 

0024 

1X0-108 

0025 

1V0-100 . 

0026 

CALL  QPLOTI 5 ) 

0027 

22-1 

0028 

OZ-IO*«SEC/3. 

0029 

IV1-80 

0030 

IX1-S0 

0031 

DO  S  I-I.4 

0032 

CALL  QPLOTI0, I X 1 • I V 1 > 

0033 

WRITE! 1 , 100 )ZZ 

0034 

I X 1  - 1 X 1 -300 

0035 

5 

ZZ-ZZ-DZ 

0036 

Ilf 

FORMAT! F6 . 1  ) 

0037 

IXI-IX0 

0038 

IV1-IV0 

0039 

DO  2  1-1,10 

0040 

CALL  QPLOT! I, IX1.IV1  ) 

0041 

IY2-IY1-600 

0042 

CALL  QPL0T!3,IX1,IV2) 

0043 

2 

IX1-IX1-100 

004  4 

IXI-IX0 

0045 

DO  3  J-1,9 

0046 

CALL  QPLOT! 1 , 1X1 , IV1 ) 

004  7 

IX2-IXI-900 

0048 

CALL  QPLOTI 3, 1X2, IV1  ) 

0049 

3 

IVI-IV1-75 

0050 

CALL  QPLOT! 6) 

0051 

1X1-0 

00S2 

1 V I -500 

0053 

DO  20  1-1,15 

0054 

IVAT-IVA! I ) 

0055 

CALL  QPLOTI 2 ,IX1,IY1,IYAT,2) 

0056 

21 

m-IYl-20 

0067 

1X1-400 

0058 

IV1-40 

0059 

CALL  QPL0T(2,IX1,IV1,IXA,14) 

0060 

1X1-50 

0061 

IV1-IV0 

0062 

DO  B  1-1,5 
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0063 

DO  7  J-l.Z 

0064 

IZ«2*( I  -  1 >♦ J 

0065 

7 

ii(  j >*m<  iz> 

0066 

call  atLom.m.m 

0067 

6 

m-ivi*is* 

0068 

CALL  QPLOT < i.IXM. 1Y» 

0069 

DO  1  1-1, ID 

0070 

IY-1Y«*1DATI I  ) 

0071 

IX-1*9Z#./JD*)M 

0072 

CALL  QPLOTI 3 , IX . IV ) 

007  3 

1 

CONTINUE 

0074 

END 

0075 

ENOS 

#CFFT  T -00003  IS  ON  CR00002  USING  MBA  SIRS  R •Ml* 
1111  FTN4.L 

0002  SUBROUTINE  CFFT(X.M) 

0003  C  WILL  OBTAIN  COMPLEX  FFT  OF  ARRAY  Y. 

0004  C  M  IS  POWER  OF  TWO  WHICH  GIVES  ORDER  OF  ARRAY 

0006  C 

0006  COMPLEX  X( 1 ),U.W,T 

0007  N-2**M 

0008  N2-N/2 

0009  N 1  - N - 1 

0010  0*1 

0011  DO  3  I-I.N1 

0012  I F  (  I  .GE.  J  >G0  TO  1 

0013  T*X( 0  > 

00U  x<o;»x<n 

0015  X< I >«T 

0016  1  K-N2 

001 7  2  IF( K  .GE.  J  >GO  TO  3 

0018  J-O-X 

0019  K-K/2 

0020  GO  TO  2 

0021  3  0 ■ J 

0022  PI-3.1415926 

0023  00  5  L- l , M 

0024  L  E -2**L 

002S  LEl-LE/2 

0026  U-( 1.0,0.01 

0027  W-CMPIXl COS (PI/LEI ),SIN(PI/LE1 ) ) 

0028  DO  5  J-l.LEl 

0029  DO  4  I-J.N.LE 

0030  I D- I >1 E 1 

0031  T-X( ID  )*U 

0032  X( ID  )*X( I  )-T 

0033  4  X(I)-X<n*T 

0034  5  U»U*W 

0035  RETURN 

0036  END 

0037  ENDS 
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#PUFT  T-BBBB3  IS  ON  CR I«ff2  USING  000V 2  BLKS  R-0000 


0001 

FTN4 , 

,  L 

0002 

SUBROUTINE  PUFT<  V ,M » I  VP ) 

0003 

COMPLEX  Y( 51 2  > 

0004 

N-2**M 

0005 

DO  1  I  - 1  ,  N 

0006 

I 

Y ( I >-<0.0, 0.0) 

0007 

IS-N/2-IVP/2M 

0008 

IH-N/2+ IWP/2 

0009 

DO  2  I-IS, 1H 

0010 

2 

Y<  1  >-< l .0,0,0) 

001 1 

CALL  CFFT(Y.M) 

0012 

RETURN 

0013 

END 

0014 

ENDS 

#CIFT 

T-00003  IS  ON  CR0000Z  USING  00003  BLKS 

0001 

FTN4 , 

,1 

0002 

SUBROUTINE  CIFT(X.M) 

0003 

COMPLEX  X< 1 ) , U ,W, T 

0004 

N-2**M 

0005 

B-l . 0/N 

0006 

DO  6  J  - 1 ,  N 

0007 

6 

X< J  >-B*X< J  > 

0008 

N2-N/2 

0009 

Nl-N-1 

0010 

0-1 

0011 

DO  3  I-1.N1 

0012 

I  F  <  I  .CE.  0  >CO  TO  1 

0013 

T-X<  J  ) 

0014 

X< J  )«X<  I  ) 

0015 

X< I >-T 

0016 

1 

K-N2 

001  7 

2 

IF  <  K  .GE.  J )GO  TO  3 

0018 

0-0-K 

0019 

K-K/2 

0020 

GO  TO  2 

0021 

3 

0  »  J  *K 

0022 

PI-3.1415926 

0023 

DO  5  L-l.M 

0024 

L  E -2**L 

0025 

LEl-LE/2 

0026 

U-< 1 .0,0.0) 

0027 

V»CMPLX<C0S(PI/LE1  ).-SIN(PI/LE!M 

0028 

DO  5  0-1, LEI 

0029 

DO  4  1-0, N, IE 

0030 

I  D-  I  ♦L  E  1 

0031 

T-X < ID  )*U 

0032 

X<  ID  ;-X<  I  )**T 

0033 

4 

X< I  )-X<  I  )*T 

0034 

5 

U-U*W 

0035 

RETURN 

0036 

END 

0037 

ENDS 

#PU  WIII3  IS  ON  CRf0»2  USING  00002  BLKS  R-0000 

0001  : TR ,  #RT 

0002  : LG . IS 

0003  : HR , XMA 1 N 

0004  : MR , XE  R  TNN 

0005  : MR , XDUM 

0006  : MR , XP  UTG 

0007  : MR , X NOYES 

0008  : MR , XP  UAQ 

0009  : MR  ,  XOF  PK 

0010  : MR , X  P  SN 

0011  : MR , XP 1  NT 

0012  MR.XPSNA 

0013  : MR . X  S  YN0 

0014  : MR . XRC VR 

0015  : MR . XF  RE  0 

0016  ; MR , XRA  T0 

0017  : MR  ,  XPHA0 

0018  : MR , XCN  T 0 

0019  : MR . X F  RO0 

0020  : MR , XOL AY 

0021  : RU.LOAOR . 99 . 1 , 18 , 1 ,2 

0022  RT  ,  LOADR 


#PUA  T-00003  IS  ON  CR00002  USING  00001  BLKS  R-0000 

0001  : TR , #RT 

0002  : l G  ,  1 0 

0003  : MR , XP  UAN 

0004  : MR , XPUANZ 

0005  : MR , XOF  UP 

0006  : MR , XP  UF  T 

0007  : MR , XC I F  T 

0008  : MR  .  XCFF T 

0009  :RU, LOADR .99, 1 , 18,0,2 

0010  :RT, LOADR 


#PUA3  T-00003  IS  ON  CR00002  USING  00001  BLKS  R-0000 

0001  : TR ,#RT 

0002  : LG , 5 

0003  : MR , XPUAN3 

0004  : MR , XQPLOT 

0005  : MR , XERAS 

0006  :RU,LOAOR.99f 1 , 18 

0007  :TR.#RT 


0020  :ST,0,XOFPK: :-3:5 
0021  :ST,8,XPSN: :-3:5 
0022  : ST , 8 , XP I  NT : : -3 : 5 

0023  :ST,8,XPSNA: :-3:5 
0024  :$T.8.XSYN0: : -3 : 5 

0025  :ST,0,XRCVR: ;-3s5 
0026  :ST,0,XFREQ: :-3:5 
0027  :STf8,XRAT0: :-3:5 
0028  :ST,8.XPHA0: :-3:5 
0029  :ST.8.XCNT0: t-3:5 
0030  :ST,8.XFRO0: :-3:5 
0031  :ST,8,%0LAY: : -3  s  S 

0032  :STt8,XPUAN: :-3:5 
0033  :ST.8.XPUAN1 : :-3:5 
0034  : ST , 8 , XOF  UP : s -3 : 5 

0035  :ST,0.XPUFT: :-3:5 
0036  tST.B.XCIFT: :-3:5 
0*J7  :ST.8.XCFFT: :-3:5 
0038  :ST, 8, XQPLOT: :-3:5 
0039  :ST,8.XERAS: :-3:5 
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mi m*i*mmm*** 


#DTOT  T  -MM3  IS  ON  CUHH  US  I  NS  MM2  81 KS  ****** 


•Ml  :ST.#PUTS,8 
•MZ  :  ST ,  *PUAQ  ,  8 
•M3  :ST,*PUAN,8 
8004  :ST,#PUANl,8 
0005  :  ST , 0PUAN2 • 8 

0006  :ST,#PUAN3,8 
0007  : ST , #PUF  T , 8 

0008  : ST , #CF  F  T  ,  8 

0009  :ST,#CIFT.8 
0010  :ST,#PU,8 

0011  :ST.*PUA.8 
00 12  : ST , 4PUA3 . 8 

00 13  :ST.«*T,0 

0014  :ST.XMAIN,8 
0015  : ST , XE  RTNN , 8 

0016  :  ST , XDUM, 8 

0017  i ST , XP  UTS , 8 
0018  : ST , XNOVE  S  , 8 

0019  : ST , XPUAQ , 8 

0020  : ST , XOFPK , 8 

0121  sST.XPSN.8 
0022  iST.XPINT.8 
••23  : ST , XP  SNA , 8 

•024  :ST.XSYN0,8 
0025  : ST , XRCVR , 0 

••26  : ST , XPRCO , 8 

0027  :ST.XRAT«,8 
0028  : ST . XPHA0 , 8 

0029  : ST , XCNT0 , 8 

0030  : ST . XFRO0 , 8 

0031  :ST,XDLAV, 8 
0032  : ST , XPUAN , 8 

0033  : ST , XPUAN1 , 8 

8834  -.ST.XOFUP.8 
0035  SST.XPUFT.,8 
••36  : ST . XC I  FT ,8 

0*37  1  ST , XCFFT , 8 

0038  : ST , XQPLOT , 8 

•039  : ST . XERAS , 8 


•TTOD  T -0000 3  IS  ON  CNMM2  USIN6  MM2  SLKS  »■»»»» 

•MI  :ST,8,#PUTC::-3 
»*»Z  : ST , 8 , #PUAO : : -3 

•••3  : ST ,8,#PUAN: : -3 

0004  tST.8.#PUANl : :-3 
•«• 5  :ST,8,*PUAN2: :-3 

0M6  :ST,B,#PUAN3:  :-3 
•007  :ST,8.#PUFT: :-3 
•008  :  ST , 8 ,#CFF  T : : -3 

•••9  :ST,8,#CIFT: ;-3 
00 1»  ! ST , 8 , #PU  J  J  *3 

0011  : ST , 8 , #P  UA : : -3 

••12  *  ST , 8 . #P  UA3 : : -3 

••13  :ST,8,#RT: :-3 

••14  :ST,8,XMA1N: :-3i5 
••15  :ST,8,XERTNN::-3iS 
••16  1ST, 8, XDUM: : -3 : S 
••17  :ST,8 ,XPUTC : : -3 : 8 

••18  1ST, 8 .XNOVE S: : -3 : 5 
••19  :ST, 8, XPUAQ: :-3:5 
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END 


DATE 

FILMED 


